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Filling transitions on rough surfaces: inadequacy of Gaussian surface models 

Renaud Dufour, Ciro Semprebon, and Stephan HerminghaufJ^] 

Max Planck Institute for Dynamics and Self-Organisation, 37077 Gottingen, Germany 

We present numerical studies of wetting on various topographic substrates, including random 
topographies. We find good agreement with recent predictions based on an analytical interface- 
displacement-type theory mm, except that we find critical end points within the physical parameter 
range. As predicted, Gaussian random surfaces are found to behave qualitatively different from non- 
Gaussian topographies. This shows that Gaussian random processes as models for rough surfaces 
must be used with great care, if at all, in the context of wetting phenomena. 


PACS numbers: 68.05.-n, 05.40.-a, 64.75.-g, 68.08.- 

I. INTRODUCTION 

While the physics of wetting and spreading on ideally 
smooth surfaces has reached a status of mature textbook 
knowledge, the wetting properties of randomly rough 
solid substrates, which are by far more relevant, are still 
poorly understood. This is in part due to the wide range 
of length scales involved. Its lower bound is the charac¬ 
teristic length scale of the molecular interactions shaping 
the three-phase contact line region, which is on the order 
of a few tens of nanometers. Topographic roughness or 
chemical inhomogeneities on length scales smaller than 
that are integrated by these forces into quantities like 
the contact angle or the contact line friction. Its upper 
bound is set by the capillary length, above which gravity 
effects become dominant. It is usually of the order of a 
millimeter (for water, it is about 2.7 mm). The lateral 
scale of typical roughness topographies of natural sur¬ 
faces is quite generally found within this range, which 
extends over five orders of magnitude. 

Despite its obvious importance, this range of random 
topographies has so far only scarcely been treated as to 
its effect on wetting. Most authors have tried to focus 
on isolated aspects of the topography using simplified 
model geometries, like rectangular [3j, triangular la , or 
algebraic grooves [S]- It has been found that at a cer¬ 
tain contact angle and vapor pressure, a filling transition 
occurs at which the formerly dry grooves are filled with 
liquid IMS]. The transition can be continuous or dis¬ 
continuous, depending on the geometry and the surface 
energy of the (regular) topography. 

Whenever the random nature of the roughness was 
taken into account, surfaces were usually modeled by 
Gaussian random processes isma, but filling transitions 
on randomly rough surfaces have been addressed only 
very recently. On the basis of an interface displacement 
model, it has been suggested that filling transitions are 
generic for a wide range of roughness occurring natu¬ 
rally mm- Corresponding phase diagrams have been 
given, based on certain approximate assumptions. Re¬ 
markably, it was predicted that it is precisely the Gaus- 
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sian surfaces which behave qualitatively different from 
all other random topographies. Since no real surface is 
precisely Gaussian, this has potentially wide ranging con¬ 
sequences. 

Given the practical importance of this prediction, 
deeper investigations are necessary to test the predictive 
power of the interface displacement theory, and thereby 
to verify its inherent assumptions. In the present paper, 
we investigate the wetting phase diagrams for a number 
of different surface topographies, including simple model 
geometries as well as random surface profiles of real 
samples. In accordance with the predictions mentioned 
above, we find that it is essentially the non-Gaussian fea¬ 
tures which gives rise to filling transitions. 

The first systematic study of wetting on an irregularly 
rough substrate at finite contact angle is due to Wenzel 
M- He characterized the roughness by a single parame¬ 
ter, r, which he defined as the ratio of the total substrate 
area divided by the projected area. The free energy which 
is gained per unit area when the rough substrate is cov¬ 
ered with a liquid is given by r( 7 Sff — "f s i ) (subscripts si 
and Ig standing for solid-liquid and solid-gas interfaces, 
respectively). If this is larger than the surface tension 
of the liquid, 7 , the macroscopic contact angle vanishes, 
because covering the substrate with liquid releases more 
energy than is required for the formation of a free liquid 
surface of the same (projected) area. More specifically, 
force balance at the three-phase contact line yields 

~ r (lsg ~ Isl) a 

cos 0 =---= r cos 6 ( 1 ) 

7 

for the macroscopic contact angle, 0 , on the rough sur¬ 
face. 9 = arccos{( 7 ss — 7 s ;)/ 7 l is the microscopic con¬ 
tact angle according to Young and Dupre. When this 
is smaller than 0w = arccos(l/r), which we will hence¬ 
forth call Wenzel’s angle, 0 vanishes, and the substrate is 
likely to be covered with an ‘infinitely’ thick liquid film. 
We focus here on the vicinity of this transition, consider¬ 
ing liquid-vapor coexistence, but also the off-coexistence 
case (i.e. below the saturated vapor pressure). 

The article is organized as follows. In section |TT] we 
briefly recall the framework of the interfacial displace¬ 
ment model on which we base our analysis. In section 
|III| we introduce the numerical method to compute equi¬ 
librium morphologies of the liquid interface on irregular 
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FIG. 1. Substrate topography g( r) and interface of the ad¬ 
sorbed liquid film /(r) (dots) of (very small) finite thickness 
lo . The microscopic contact angle is indicated by 6. 


surfaces. In section |IV[ we consider a simple periodic 
surface topography and compare the numerical results 
to the results of the interface displacement theory. In 
section [V] both methods are applied to more complex, 
irregular substrates in order to confirm the generic char¬ 
acter of the filling transition. Finally, in section [VTl we 
apply the interfacial displacement model to investigate 
the filling transition of a real physical roughness, as well 
as one having nearly Gaussian properties. 


II. INTERFACE DISPLACEMENT MODEL 

Let us briefly review the interface displacement the¬ 
ory which has been put forward before (2j. The aim of 
the theory is to predict the amount of adsorbed liquid 
as function the microscopic contact angle and the mean 
curvature of the liquid interface. 

We consider a homogeneous and isotropically rough 
solid surface defined on a set f 1 in the plane by a func¬ 
tion z — f( r), where z is the height of the surface above 
the point r = ( x,y ). We assume that V/ and A/ exist 
everywhere, and that the roughness amplitude is much 
smaller than its lateral length scale, V/ < 1, which is 
the case for a large class of natural roughness. The free 
surface of a liquid film adsorbed on this substrate is de¬ 
scribed by a second function, g( r), which deviates ap¬ 
preciably (i.e., by more than Iq) only in the wet area, 
We SI. The morphology of the liquid/vapor interface, 
and thereby the amount of adsorbed liquid, is determined 
by the boundary conditions 

g( r) = /(r) + l 0 Vr € dW 
|V( 5 (r)-/(r))|«0 VredW (2) 

where 5 - / > lo on the wetted patches of the sample. 
The last equation is to be fulfilled to first order in 0, V/ 
and V g everywhere on <9W\ which is the projection of 
the three-phase contact line. Applying Green’s theorem 
to (<7 — /), we can write 

/ n • V(t/ — /) ds = [ A(g — /) d 2 r (3) 
JdW JW 

where s is the distance along dW and n its unit normal 
vector. Writing Eq. ([2| as n ■ V (g — /) ss 6 on dW and 


combining with Eq. (|3|, we obtain 

W + [ [2H — A/] d 2 r = 0 (4) 

Jw 

where l denotes the length of dW per unit surface area 
and H is the mean curvature of the liquid/vapor inter¬ 
face. In order to evaluate Eq. Q, it is worth describ¬ 
ing a given topography in terms of statistical quantities. 
Therefore we define: 


p(z) = (fll - 1 / n S(z - /(r)) d 2 r (5) 

°ioo =p(^)" 1 i^r 1 /ni v /( r )i <5 ( 2; -/( r )) d2r ( 6 ) 

a 2 (z) = p(~) _ 1 |fl |~ 1 J n |V/(r )| 2 S(z — /(r)) d 2 r (7) 


where p(z) is the surface height distribution, cri(z) and 
a 2 (z) describe the average slope (|V/|)^ and average 
square slope (|V/| 2 ) Z , respectively, at elevation 2 . Equa¬ 
tion Q can then be re-written as [ 2 ]: 


2 H f p(f)df 

J —OO 


~ 0~2(z) 

.^lO) 


°i{z)p{z) (8) 


We finally define the cumulative height distribution, 
W(z) = J^ oo p(f)dz and note that the length of a con¬ 
tour line at elevation h is given by L(h) = cn f h Wh)[T5l. 
It is one of the major approximations used previously 
mm to assume that L(h) is roughly equal to the length 
of the contact line if h is taken to be the average height 
of that line, and that W ( h ) is approximately equal to the 
wetted sample area. We then obtain 


2 HW{h) 


~o- 2 (fr) 
.oi (h) 


L(h). 


(9) 


If p(z), <J\{z) and a 2 (z) are known from experimental 
characterization of the sample, Eq. ([9]) allows to de¬ 
termine the equilibrium level h of the liquid filling the 
troughs, as a function of the microscopic angle 6 and 
curvature H. A reasonable estimate of the volume of 
adsorbed liquid per unit area can then be obtained by 

V(h)= ( ( h~z)p(z)dz (10) 

J —OO 


note that such approximation neglects the curvature of 
the liquid meniscus with respect to the horizontal plane 
defined by the cut at the level 2 . Nevertheless, the be¬ 
haviour of V thus defined will qualitatively reflect all fea¬ 
tures of interest of an adsorption isotherm. 


III. NUMERICAL ENERGY MINIMIZATIONS 

In order to assess the accuracy of the interface dis¬ 
placement model, we computed equilibrium profiles g( r) 
of the free liquid interface from numerical minimizations 
of the free energy: 

H9{r)}=\ [ d 2 r |V<?(r )| 2 

J 

+ f d 2 r <I> (£(r)) — A f d 2 r £(r) , 

Jn Jn 


( 11 ) 
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with the interfacial tension of the liquid-air interface, 7 , 
the local thickness of the liquid film, £(r) = g( r) — /(r), 
and the domain of integration, fi. A sketch of a cross 
section perpendicular to the substrate is shown in Fig. |Tj 
The parameter A in front of the last term on the RHS 
of Eq. © represents the Laplace pressure of a reservoir 
that can exchange liquid with the substrate, i.e., A = 
2 ~fH. 

To obtain equilibrium profiles with a desired appar¬ 
ent contact angle, we use a generic short ranged effective 
interface potential, 


$(0 = $0 


24 


£(1 + VS) 


- 1 exp - 


£(1 +V5)' 


24 


( 12 ) 


The equilibrium film thickness 4 corresponds to <f>'(Zo) = 
0 (vanishing disjoining pressure). The depth of the global 
minimum of interface potential determines the apparent 
contact angle cos 9 = 1 — | $( 4 ) | h- 

For a given substrate topography /(r), the free energy 
(Eq. (Ill) is numerically minimized employing the public 
domain software Surface Evolver ehj. To this end, the 
free interface is represented as a mesh of small triangles 
the nodes of which are moved with a conjugate gradi¬ 
ent algorithm. The triangulation is regularly re-meshed, 
such as to keep the typical size of the edges below about 
~ A/20, where A is the shortest wavelength of the to¬ 
pography, and to avoid the formation of long and narrow 
triangles. 4 was set equal to A/20. It was found that the 
numerical stability could be further improved by intro¬ 
ducing an additional hard repulsive constraint ensuring 

£ > 4- 

In a typical run, the computation is initialised by 
bringing the liquid interface in touch with the highest 
peak of the roughness. The contact angle 0 is maintained 
constant, while a small but finite pressure is applied. By 
ramping up the pressure parameter A (Eq. ©) , we ex¬ 
plore sequences of interfacial configurations of the liquid 
with a continuously varying mean curvature, H. The 
physical reason for this preference is that if one starts 
instead from a small amount of adsorbed liquid, core¬ 
sponding to large interfacial curvature, the liquid consists 
of a large number of isolated liquid puddles, and it is not 
straightforward how to generate or select a well defined 
initial state. 

Note that the sequences of morphologies obtained with 
this procedure do not necessarily represent global energy 
minima. This is resembles an experiment where the par¬ 
tial pressure of the vapour of the adsorbed liquid is grad¬ 
ually reduced. Given the fact that the typical length 
scales of natural roughness correspond to interfacial en¬ 
ergy barriers much larger than thermal energies, global 
free energy minima are reached in experiments only af¬ 
ter close to infinite times. It should be noted that a 
similar approach was already successfully adopted to in¬ 
vestigate the instabilities of a liquid meniscus leading to 
advancement of the liquid front on a regularly patterned 

substrate [13- 



fig. 2. Statistical properties of the regular topography from 
Eq. ( |13| ) . From top to bottom: (j) = (— n/2,(f) = 0,7r/2). 
Solid curves represent the height distribution p(z). Condi¬ 
tional mean slope cr 1 (z) and a 2 (z) are shown as dotted and 
dashed lines, respectively. Inset represents the unit cell of the 
periodic surface. 


IV. EGG-CARTON TOPOGRAPHY 


The major difference between one dimensional periodic 
topographies and truly two dimensional substrates is the 
presence of saddle points in between peaks and valleys in 
the latter. To assess the fundamental mechanisms under¬ 
lying the drying and filling transition on two-dimensional 
topographies, we begin our analysis considering a sim¬ 
ple example. It consists of a hexagonal arrangement of 
peaks, valleys and saddle points, whose locations can be 
set by tuning appropriate parameters. We will refer in 
the following to it as ‘egg carton’. An hexagonal egg- 
carton can be mathematically described by 

( . 4ttui . 47 tu 2 . . Anu 3 \ 

z = Al sin(—-1- (j)) + sm(——) + sm (——) 1 (13) 

where L is the unit length of the periodic cell, 
A = 0.03 x L tunes the amplitude of the corrugation 
and Ui represent hexagonal coordinates, 


u 1 


u 2 


u 3 


= x 


x + V3y 
2 

x - V3y 


(14) 


2 
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FIG. 3. Graphical solution of Eq. § for the egg-carton pat¬ 
tern with (j> = 0 (Fig. [5] middle panel) and two combinations 
of [9, H] values. Left : 9 = 5° and H — 1. This configura¬ 
tion admits a single solution at z ~ 0.025. Right : 9 = 20° 
and H = 0.4. This configuration admits multiple solutions 
of which two are stable (closed circles) and two are unstable 
(open circle). The stability of the solution depends upon the 
slopes of the intersecting curves. 


In Eq. (13), the phase variable cj) enables to tune the verti¬ 
cal level of the saddle point. We consider three different 
configurations, with <fi = —7t/2, (f> = 0, and 4> = tt/ 2. 
These choices place the saddle point close to the top, in 
the mid-plane, and close to the bottom of the topography, 
respectively. It can be observed in Fig. [2] that the level of 
an isolated saddle point significantly affects the shape of 
the distributions p(z), (J\ (z) and cr 2 (z) (cf. Eq. 0 , @). 
and 0, respectively). 

The average level h of the contact line of an adsorbed 
liquid film can be predicted with the interface displace¬ 
ment approach by solving Eq. ([ 9 ]). For given values of 
9 and the mean curvature of the free liquid surface, H, 
the solution can be obtained graphically. To illustrate the 
mechanism employed in the numerical solution of Eq. ([ 9 ]) , 
let us consider the case </> = 0 in the examples depicted 
in Fig. [3] Depending on 9 and H, Eq. 0 admits either 
single or multiple solutions. The stability of the solu¬ 
tions given by the graphical intersection can be inferred 
by considering that Eq. 0 represents a force balance. 
For a solution to be stable, upon a displacement of the 
three-phase contact line, the unbalance of wetting forces 
have to bring the film back toward the equilibrium point. 


Fig.[4ji shows a series of adsorption isotherms, in terms 
of the average level of the contact line h as function of 
the curvature H. They demonstrate satisfactory agree¬ 
ment between the interface displacement approach (solid 
curves) and the numerical energy minimizations (dashed 
curves) for all cases considered. For a saddle point lo¬ 
cated either close to the peaks or to the valleys of the 
topography (<j) = ±7t/ 2), the desorption curves exhibit a 
single discontinuous jump from a solution of finite film 
thickness to the dry state. This transition reflects the 
fact that the topography cannot accommodate the pre¬ 
vious liquid configuration anymore at the given H and 9 
on the substrate topography. The transition line is anal¬ 
ogous to the off-coexistence filling transitions predicted 
for simple translationally invariant substrates mm- 
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FIG. 4. a) Comparison of the desorption curves obtained for 
the egg-carton pattern from the theory (continuous lines, sta¬ 
ble solutions only) and energy minimization (dashed lines) for 
the three morphologies shown in Fig. [5] h* is the film level 
normalized w.r.t. the rms roughness, and H* is the film cur¬ 
vature normalized w.r.t. the rms curvature. Grey areas rep¬ 
resent values outside the upper (h+) and lower (/i_) bounds 
of the patterns. Results are shown for 9 = 5° (black), 15° 
(dark grey), and 25° (light grey), b) Wetting phase diagrams 
(same parameters as in the three top panels). The transition 
lines from the theory and energy minimizations are shown as 
continuous lines and squares, respectively (dashed lines are 
guides to the eye). The transition lines are obtained by lo¬ 
cating the discontinuities in the desorption curves, and end 
at critical end points, as indicated by the closed symbols. 


The largest allowed curvature corresponds to the cur¬ 
vature of the topography in the single periodic valley of 
the egg-carton, and is attained for vanishing 9. In this 
case the adsorption isotherm exhibits a continuous be¬ 
havior. Note that in the case <j> = — 7t/ 2, the liquid film 
consists of a single connected domain, while in the case 
(j) = 7t/ 2, the liquid film consists of a periodic arrange¬ 
ment of identical isolated domains. As long as as the 
statistical properties are similar, our analysis shows that 
for a periodic substrate the topology of the liquid do¬ 
main plays no role in the adsorption or desorption. The 
case 4> = 0, where the saddle point is located around 
z = 0, introduces a relevant qualitative difference, as 
drying proceeds with a change of topology when the liq- 
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uid interface dewets the saddle point. Consequently, we 
observe an additional discontinuity, splitting the desorp¬ 
tion curve into two branches, both characterized by a 
finite amount of adsorbed liquid. Note that this feature 
is captured in both the interface displacement approach 
and the numerical calculations. 

It is worth noting that the interface displacement ap¬ 
proach captures the presence of the two branches, be¬ 
cause the statistical properties are strongly affected by 
the presence of the saddle point (cf. Fig. [2]). On the 
other hand, a close look at the numerically computed 
shape of the liquid interface shows that the contact line 
experiences a stronger deviation in the vertical direction 
in the vicinity of the saddle point. Such configuration de¬ 
parts from the approximation that length of the contact 
line is approximated by the length of surface isolines, as 
required by the interface displacement approach to ob¬ 
tain Eq. ([ 8 ]) . Consequently, the agreement is slightly less 
accurate. 

Fig# shows the resulting phase diagrams in the plane 
spanned by 9 and H. The solid curves are obtained from 
the interface displacement model, while the squares re¬ 
present the discontinuities in the curves displayed in the 
top panel. Obviously, the qualitative agreement is very 
good, as the number of observed discontinuities is the 
same in all three bottom panels. The quantitative agree¬ 
ment concerning the intercept and slope of the phase 
boundaries is satisfactory as well, given the crudeness 
of the approximations used in the interface displacement 
model. 

So far we can conclude that the interface displacement 
model put forward earlier m is capable of reproducing 
the structure of the adsorption isotherms rather well, al¬ 
though substantial simplifications have been involved in 
the derivation of its central equation 


V. IRREGULAR ROUGHNESS 

Let us now turn to more irregular roughness topogra¬ 
phies. For the numerical investigations, we have to ensure 
that the boundary conditions are met at the limits of the 
support. This can be best achieved by imposing periodic 
boundary conditions. Consequently, we generated peri¬ 
odic roughness profiles /(r) by means of a superposition 
of many (> 400) sinusoidal modes shifted by randomly 
chosen phases 0 k £ { 0 , 27r} 

/(r) = ^ a k sin (k • r + 0 k ) (15) 

k 

with r = (cc, y) £ U = [ 0 , 1] x [ 0 , 1] and two dimensional 
wave vectors k = 27t where m,n £ {— 10 ,..., 10 }. 
The real Fourier amplitudes decay according to a power 
law, 


n k 


0 for k = 0 

A 0 (l + |k 2 |) -1 / 2 else 



FIG. 5. Example of desorption curves obtained from the 
numerical energy minimization for e = 0 (a), e = 0.02 (b) 
and e = 0.05 (c). V* is the adsorbed volume per unit area, 
normalized with the rms roughness amplitude. In each panel, 
the contact angle starts at 91° and increases in one-degree 
steps. The black portions of the curves connect the desorption 
branches through the discontinuity (i.e. circles reported in 
Fig# 


This leads to topographies with nearly Gaussian proper¬ 
ties (see below). 

In order to create topographies with non-Gaussian 
properties, we distort the so-obtained profile by apply¬ 
ing the transformation 

/(r) = /(r) + e-/(r) 2 . (16) 


The effect of Eq. (16) is to sharpen the peaks and flatten 
the valleys of / as compared to /. For e = 0 the pattern 
is unchanged and has (nearly) Gaussian properties. For 
e > 0 we depart from Gaussian properties as the height 
distribution becomes skewed, and a positive correlation 
is introduced between the elevation 2 and surface slope. 
A measure for the deviation of the chosen profile from 
Gaussian properties is the third moment (skewness) Sk 
of the height distribution, which vanishes for a perfectly 
Gaussian topography (20] [2T]. We studied the cases e = 


0, e = 0.02, and e = 0.05. 

A typical set of results for the energy minimizations 
is reported in Fig. [5] In the case of a nearly Gaussian 
topography (e = 0 ) we observe only a few discontinu¬ 
ities for a finite curvature, and only for an intermediate 
range of material contact angles. These jumps involve 
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FIG. 6. (a) Height distribution for the numerically generated random roughness with t from 0 to 0.05 (cf. Eq. fief), (b-d) 
Comparison between the theory and numerical results. We plot the phase diagram spanned by the normalized curvature H* 
and the ratio 9/9w, with 6w = 5.2°, 5.3° and 5.4° for panels b-d respectively. From the numerical results, the diameter of 
each circle is proportional to the height of the jump. The open circle, for scale, corresponds to a 100% jump. Note that only 
events assessing for at least 10% of the total volume are shown. The continuous black line is the transition predicted by the 
interface displacement model. 


also a relatively small volume variation, indicating that 
after the transition many troughs remain filled with liq¬ 
uid, and continuously dry out until the largest curvature 
allowed is reached. This finding can be related to the 
presence of saddle points at any level of the topography, 
whose cumulative effect involves a dramatic change in 
the topology of the liquid interface from a single or a 
few large domains to a large number of small domains. 
Both for larger and smaller contact angle we observe a 
continuous drying, but while for small 9 the drying starts 
from a completely filled roughness, for large 9 the initial 
state, for zero curvature, involves already a partially dry 
topography. 

On the other hand, for the case e = 0.05, we observe a 
large number of jumps at intermediate curvatures, most 
of which lead directly to the dry solution. This find¬ 
ing shows that the stability of the upper branches of the 
liquid morphologies is reduced when most of the saddle 
points are located close to the bottom. 

In contrast to the egg-carton topography, the location 
of discontinuities from numerical energy minimizations 
on an irregular topography, as in the present case, can 
be affected by strong fluctuations related to local details. 
To improve comparability, we therefore created ten sam¬ 
ple topographies by assigning different series of random 


phases, and performed ensemble averaging. The inter¬ 
face displacement model is then applied to the equiva¬ 
lent pattern obtained by combining the statistical prop¬ 
erties of all ten samples. For e = [0,0.02,0.05], the 
corresponding ensemble averages for the skewness are 
Sk = [0.006,0.330,0.790]. The results from the energy 
minimizations instead are analysed by detecting the cur¬ 
vature H and the magnitude AV of the volume differ¬ 
ence. Each jump is then represented by a circle in Fig. [ 6 j 
centered at the point where the jump was found, and 
with a diameter proportional to AV. Note that the con¬ 
tact angle is normalized by the respective Wenzel angle 
9w for each topography. 

As we can observe in Fig. [Gj the position of the tran¬ 
sition line predicted by the interface displacement ap¬ 
proach is consistent with the discontinuities identified 
with the energy minimizations. For the case e = 0 we 
observe a short flat line, surrounded by a few spots of 
small radius. Contrarily for the case e = 0.05 we ob¬ 
serve a clear transition line, surrounded by large spots. 
Remarkably, Fig. 0' -d shows that in the limit H* -4 0, 
the transition lines and numerical events depart from a 
contact angle 9f larger than the Wenzel angle 9w- This 
behavior is in agreement with the interface displacement 
theory [2] (see also discussion below). In particular, for 
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the nearly Gaussian configuration (e = 0), we expect 
df/dw ~ Z/Vn ~ 1-13, in very good agreement with 
Fig. [BJd- 


VI. NEARLY GAUSSIAN VS REAL 
ROUGHNESS 


In the previous section, we validated the interface dis¬ 
placement model for a complex irregular roughness, but 
we were still limited to a relatively small number of cor¬ 
rugation modes. In the present section, we apply the 
interface displacement approach to more complex sub¬ 
strates, in order to compare a nearly Gaussian roughness 
to the real topography of a sand blasted surface. 

Let us first consider an ideal Gaussian roughness. We 
assume the profile to have vanishing mean, rms rough¬ 
ness a, and correlation length r. The height distribution 
is then given by p(z ) = Af(0 , a 2 ), where A7(0, a 2 ) denotes 
the Gaussian normal distribution with variance a 2 and 
zero mean. Additionally, ay = ^/na/r and ay = 4(a/r) 2 
are constants mm, indicating that the probability of 
finding a certain slope at a given level, z, is independent 
of z. This is linked to the fact that an ideal Gaussian 
roughness has an infinite support. Consequently, Eq. (19]) 
admits a unique solution for 9 < Of = ^ ( \/H ), 

and no solution for 9 > 9f. Therefore the phase dia¬ 
gram of a ideal Gaussian roughness is characterized by a 
horizontal line at 9 = 9f delimiting the transition from 
wet to dry configurations. It is instructive to compare 
9f with the Wenzel angle 9\y- For Gaussian roughness, 
the Wenzel parameter is r « 1 + (a/r) 2 [22] , This leads 
to 9\y ss 2 a/r and to a fixed ratio 9f/9w ~ 2/-^. For 
a non-Gaussian roughness, we also expect 9f to be sys¬ 
tematically larger than 9w [2], as we already observed 
in the previous section (cf. Fig. [Gjo-d) . In practice, this 
means that if the liquid does not wet the substrate well 
enough to fulfill the Wenzel condition for complete wet¬ 
ting (Eq. 0). it may nevertheless form a wetting layer. 
In the present case, we expect the transition from dry 
to wet being dominated by complex energy barriers as 
described for wicking and imbibition phenomena |2~5] . 

We investigate the numerical solution of Eq. ([9]) on 
nearly Gaussian roughness, by limiting the amplitude of 
the topography in the real space. The surface is obtained 
by taking the inverse Fourier transform of an isotropic 
Gaussian power spectrum, 


S(k) 


2 


|k| exp 




(17) 


where k = (k x , k y ) is the 2D wave vector. The rms am¬ 
plitude a and correlation length r are set to 2.3 /mi and 
21.5 /tm, to match the properties of the topography of 
the the experimental sample. The latter is obtained by 
sandblasting of a copper plate, with sand grains of size 
in the range 70 — 110/tm (Rohde AG). The surface pro¬ 
file is measured with a white light interferometer (Weeko 


NT1100) over an area of 450*600 /im, with lateral and 
vertical resolutions of Sx « 0 . 8 /tm and Sz < 1 nm re¬ 
spectively. The map of the topography is pre-processed 
by applying a 5 * 5 pixels Gaussian filter in order to re¬ 
move measurement artifacts. The statistical properties 
are computed from profiles obtained at three different 
positions on the sample. 

Both the surface profiles and their statistical proper¬ 
ties p{z) and cr 12 (z) are shown in Fig. l] panels a) and 
b). As expected for the Gaussian roughness, the height 
distribution matches the Gaussian distribution, and 01,2 
appear constant over most part of the range of z , in agree¬ 
ment with expected values of ay = 0.190 and a 2 = 0.046 
respectively. Because of the limited support, there is 
a highest elevation at z = z+ and a deepest trough at 
z = z_. As a consequence, both oy and oy drop to zero 
as z -A z + /_. The height distribution of the real sam¬ 
ple seems at first glance close to the analytical Gaussian 
curve. However, a difference is found in the value of 
the skewness, Sk ss 0.15, indicating a predominance of 
valleys. This departure from Gaussian properties is also 
reflected in the positive correlation between the slope and 
the elevation. Such feature can be expected if one consid¬ 
ers the roughening mechanism underlying sandblasting, 
which is known to produce positively skewed topogra¬ 
phies [24H26] . The desorption curves obtained from the 
numerical solution of Eq. are shown in Fig.[7]panels c) 
and d). The adsorbed volume per unit area is normalized 
by the rms amplitude a. 

Two noticeable differences arise between a nearly 
Gaussian roughness on finite support and an ideal one. 
For a nearly Gaussian roughness we observe that the 
adsorbed volume at coexistence is finite V*{H = 0) = 
V*{z + ) and the topography completely dries for a fi¬ 
nite value of H*. Instead, for an ideal Gaussian one 
we would have V*(H = 0) -a 00 , while the complete 
drying would occur for H* -A 00 , because of the non 
vanishing probability of finding troughs of any size and 
curvature. Regardless of these two aspects, the other 
features are present. Fig. [7]: exhibits no bifurcation and 
the desorption curves continuously decrease until reach¬ 
ing the curvature of the smallest trough. These results, 
compared to the nearly Gaussian topography considered 
in the previous section (Fig. suggest that the few dis¬ 
continuities observed were an effect of the limited number 
of modes considered m- Increasing the contact angle, 
the largest value of 9 for which Eq. Q admits a solu¬ 
tion is 13°. This corresponds to the innermost curve in 
Fig#. If 9 is slightly increased to just 14°, no adsorp¬ 
tion is observed anymore. Hence we can identify a critical 
value of 9 above which the adsorbed amount of liquid goes 
abruptly to zero, irrespective of the mean curvature, H. 
In the case of Gaussian roughness this value has an ana¬ 
lytical expression. Its value is Of = = ^ 7=7 = 13.8°. 

The predictions of the interface displacement model 
applied to the sandblasted topography are shown in 
Fig#. Consistently with the measured statistical prop¬ 
erties of the topography, the desorption curves behave in 







FIG. 7. Statistical properties of a random roughness with Gaussian power spectrum (a) and of a real rough sample treated 
by sandblasting (b). Height distribution p(z) are shown as open circles and overlays the analytical Gaussian curve in black. 
Conditional mean slope 0 - 1 ( 2 ) and 0 - 2 ( 2 ) are shown as dotted and dashed lines respectively. Inset represents the unit periodic 
cell in (a) and a 400 x 500 pm window in (b). (c) and (d) show the corresponding desorption isotherms. Continuous and 
dotted curves correspond to stable and unstable solutions to Eq. respectively. Insets are example of graphical resolution 
with closed/open circles corresponding to stable/unstable solutions and mapped onto the corresponding curves. 


agreement with the findings presented in sec. VI for the 
case e = 0.05. For small values of 0, the model predicts 
a continuous decrease of V with H*. As 8 increases, ap¬ 
proaching the critical value 9 C « 15°, a bifurcation point 
emerges at H* « 0.06. For H* > H* no solution is 
detected, while for H* < H* both stable and unstable 
branches of Eq. ®> appear. The largest material angle 
admitting a stable solution in the case of the real rough¬ 
ness is Of ss 21° (inner left curve in Fig. 0 1 )' Finally, we 
present in Fig. [8]the phase diagram obtained for the cases 
displayed in Fig. [TJ summarizing all features described 
above. Remarkably, for the real roughness we find that 
the critical end point [ H *, 0 C ] lies in the physical region 
of the parameter space. 


ically, we corroborated that filling transitions are generic 
on random roughness topographies, and the phase dia¬ 
grams are in reasonable quantitative agreement with the 
predictions of the analytic theory. Furthermore, our re¬ 
sults support the startling prediction that Gaussian ran¬ 
dom topographies behave qualitatively differently. First, 
we find no filling transitions in adsorption isotherms on 
perfectly Gaussian roughness. Second, there is a filling 
transition when 8 is varied, but its position is indepen¬ 
dent of H. This confirms that modeling topographic 
roughness with a Gaussian random process may miss im¬ 
portant physical aspects of wetting. 


VII. CONCLUSIONS 

In this paper we reported on numerical investigations 
of wetting on various roughness topographies. Our re¬ 
sults compare favourably with predictions made earlier a 
purely analytic interface displacement approach. Specif¬ 
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FIG. 8. Phase diagram for the cases displayed in Fig. [7] 
The horizontal dashed line indicates the critical angle 9 C = 
13.5(±0.5)° which can be deduced from for the Gaussian pat¬ 
tern (Fig. [7^ and[7J;). The continuous black line describes the 
discontinuous transition as introduced in Ref. [l], obtained 
for the sandblasted surface ((Fig. and[7|i)). Crossing this 
line leads to a discontinuous jump in the amount of adsorbed 
liquid. 


[1] S. Herminghaus, Eur. Phys. J. E 35 , 43 (2012) 

[2] S. Herminghaus, Phys. Rev. Lett. 109 , 236102 (2012) 

[3] R. Seemann, M. Brinkmann, E. J. Kramer, F. F. Lange, 
and R. Lipowsky, Proc. Nat. Acad. Sci. USA 102, 1848 
(2005) 

[4] K. Rejmer, S. Dietrich, and M. Napiorkowski. Phys. Rev. 
E 60 , 4027 (1999) 

[5] C. Rascon and A. O. Parry, Nature 407 , 986 (2000) 

[6] A. Malijevsky and A. O. Parry, Phys. Rev. Lett. 110 , 
166101 (2013) 

[7] A. O. Parry, A. Malijevsky, and C. Rascon, Phys. Rev. 
Lett. 113 , 146101 (2014) 

[8] A. Malijevsky and A. O. Parry, Phys. Rev. E 91 , 052401 
(2015). 

[9] B. N. J. Persson, O. Albohr, U. Tartaglino, A. I. Volok¬ 
itin, and E. Tosatti, J. Phys. Condens. Matter 17 , R1 
(2005) 

[10] N. Savva, S. Kalliadasis, and G. A. Pavliotis, Phys. Rev. 
Lett. 104 , 084501 (2010) 

[11] D. Chakraborty, N. N. Dingari, and S. Chakraborty, 
Lan g muir 28 , 16701 (2012) 

[12] R. David and A. W. Neumann, Langmuir 29 , 4551 
(2013) 

[13] F. Bottiglione and G. Carbone, J. Phys. Condens. Matter 
27,015009 (2015) 

[14] R. N. Wenzel, Ind. Eng. Chem. 28, 988 (1936) 

[15] M. S. Longuet-Higgins, Phys. Trans. R. Soc. A 250, 157 
(1957). 

[16] K. A. Brakke, Phil. Trans. R. Soc. A 354, 2143 (1996) 


[17] C. Semprebon, S. Herminghaus, and M. Brinkmann, Soft 
Matter 8, 6301 (2012) 

[18] C. Rascon, a. O. Parry, and a. Sartori, Phys. Rev. E 59 , 
5697 (1999), 9902070' 

[19] A. Rodriguez-Rivas, J. Galvan, and J. M. Romero- 
Enrique, J. Phys. Condens. Matter 27 , 035101 (2014) 
1408.1013 

[20] W. H. Press, S. A. Teukolsly, W. T. Vetterling, and 
B. P. Flannery, Numerical Recipes: The Art of Scientific 
Computing, 3rd ed. (Cambridge University Press, 2007) 
Chap. 14, p. 994. 

[21] N. Tayebi and A. A. Polycarpou, Tribol. Int. 37 , 491 
(2004) 

[22] M. Nosonovsky and B. Bhushan, J. Phys. Condens. Mat¬ 
ter 20 , 225009 (2008) 

[23] C. Semprebon, P. Forsberg, C. Priest, and 
M. Brinkmann, Soft Matter 10 , 5739 (2014) 

[24] D. J. Whitehouse, Characterization of Solid Surfaces. 
Chapter 3, edited by P. F. Kane and G. B. Larrabee 
(Springer US, Boston, MA, 1974). 

[251 B. Griffiths, Enqineerinq Drawinq for Manufacture (El¬ 
sevier, 2003) pp. 111-133. 

[26] T. Monetta and F. Bellucci, OJRM 01 , 41 (2012) 

[27] For the patterns investigated in section [V] statistics on 
the ensemble for e = 0 agree with Gaussian proper¬ 
ties. However because of the small sample size, individual 
samples can still have have non-negligible skewness and 
exhibit jumps in the desorption curves. This is for ex¬ 
ample the case in Fig. [5ji where for e = 0 the sample 
skewness is 0.19. 







